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Abstract 

Several nontrivial properties are shown for the mean square radius of gyration R 2 K of ring poly- 
mers with a fixed knot type K. Through computer simulation, we discuss both finite-size and 
asymptotic behaviors of the gyration radius under the topological constraint for self-avoiding poly- 
gons consisting of N cylindrical segments with radius r. We find that the average size of ring 
polymers with a knot K can be much larger than that of no topological constraint. The effective 
expansion due to the topological constraint depends strongly on the parameter r which is related 
to the excluded volume. The topological expansion is particularly significant for the small r case, 
where the simulation result is associated with that of random polygons with the knot K. 
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I. INTRODUCTION 



A ring polymer is one of the simplest systems that have the effect of topological entan- 
glement. The topological state of a ring polymer is given by a knot, and it is fixed after the 
ring polymer is formed. The entropy of the ring polymer with the fixed knot is much smaller 
than that of no topological constraint. Thus, there should be several nontrivial properties in 
statistical mechanics of ring polymers with a fixed topology. Furthermore, some dynamical 
or thermodynamical properties of ring polymers under topological constraints could also be 
nontrivial. In fact, various computer simulations of ring polymers with fixed topology were 
performed by several groups |TJ, El ^|, |], FjL |5], [7], R, [TT], O, |T^]. However, there are still many 
unsolved problems related to the topological effect, such as the average size of a knotted 
ring polymer in solution. 

In the paper, we discuss how the excluded volume controls the topological effect on the 
average size of ring polymers in good solution. As a model of ring polymers we employ a 
model of self-avoiding polygons consisting of cylindrical segments with radius r. Through 
numerical simulation, we investigate the mean square radius of gyration of cylindrical self- 
avoiding polygons with radius r [|1^, |14], H|. By changing the cylinder radius r, we modify 
the excluded volume effect. Thus, we can investigate the topological effect systematically 
through the simulation of cylindrical self-avoiding polygons for various values of cylinder 
radius r. 

Let us consider the two cases when the radius r is very large or very small. When the 
radius r is very large, the simulation should be related to that of the self-avoiding polygons 



on the lattice (TT|, On the other hand, when the radius r is very small, it is related 
to random polygons with a fixed topology, as we shall see explicitly through the data. In 



fact, there is quite an interesting suggestion |16| [L7|, [18] that under a topological constraint 
the average size of ring polymers with no excluded volume should be similar as that of ring 
polymers with the excluded volume, since nontrivial entropic repulsion should be derived 
from the topological constraint. According to the suggestion, the average size of random 
polygons with the trivial knot should be given by N USAW with respect to the number N of 
polygonal nodes, where vsaw is the exponent of self-avoiding walks. Thus, the small r case 
of the simulation in the paper should be important also in the study of the topological effect 
on random polygons. 
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The outline of the paper is given in the following. In Sec. II we explain self-avoiding 
polygons (SAPs) consisting of cylinder segments. We also discuss the effective exponent 
of the mean square radius of gyration under no topological constraint R 2 . In Sec. Ill, we 
discuss various nontrivial finite-size properties of the mean square radius of gyration R 2 K for 
cylindrical SAPs with a given knot type K. The ratio R 2 K / R 2 expresses the effective expan- 
sion due to the topological constraint. Through the simulation, we find that the topological 
effect is important particularly in the small r case for cylindrical SAPs. Furthermore, the 
effective topological expansion is controlled by the parameter r. In Sec. IV, we discuss the 
asymptotic expansion of the ratio R 2 K j R 2 with respect to the number N. Finally, in Sec. V, 
we graphically explain the effective expansion of the cylindrical SAPs under the topological 
constraint, through the graphs in the N — r plane. 



II. CYLINDRICAL SELF- AVOIDING POLYGONS 

A. Cylindrical ring-dimerization algorithm and random knots 

Let us introduce a model of ring polymers in good solution. We consider self-avoiding 
polygons consisting of N rigid impenetrable cylinders of unit length and diameter r: there is 
no overlap allowed for any non-adjacent pairs of cylindrical segments, while next-neighboring 
cylinders may overlap each other. We call them cylindrical self-avoiding polygons or cylin- 
drical SAPs, for short. The cylinder radius r can be related to the stiffness of some stiff 
polymers such as DNAs || |14| . 

In the simulations of the paper, we have constructed a large number of cylindrical SAPs 
by the cylindrical ring-dimerization method |T^|. The method is based on the algorithm of 
ring-dimerization [|J, and very useful for generating long self-avoiding polygons (for details, 
see Ref. jnj). Here we note that another algorithm is discussed in Ref. || for the model of 
cylindrical SAPs, where self-avoiding polygons of impenetrable cylinders with N < 100 are 
constructed in association with knotted DNAs E0J. 

In the cylindrical ring-dimerization method, a statistical weight is given to any self- 
avoiding polygon successfully concatenated. Thus, when we evaluate some quantity, we 
take the weighted average of it with respect to the statistical weight. Some details on the 



statistical weight of successful concatenation is given in Ref. |14j| . Hereafter in the paper, 
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however, we do not express the statistical weight, for simplicity. 

Let us describe the processes of our numerical experiments. First, we construct M samples 
of cylindrical SAPs with N nodes by the cylindrical ring-dimerization method. We put 
M = 10 4 . Here we note that various knot types are included in the M random samples. 
Second, we make knot diagrams for the three-dimensional configurations of cylindrical SAPs, 
by projecting them onto a plain. Then, we calculate two knot invariants A^(— 1) and V2(K) 
for the knot diagrams. Third, we select only such polygons that have the same set of values 
of the two knot invariants, and then evaluate physical quantities such as mean-squared 
gyration radius for the selected cylindrical SAPs. 

The symbol A^(- 1) denotes the determinant of a knot K, which is given by the Alexan- 
der polynomial A(t) evaluated at t = —1. The symbol V2(K) is the Vassiliev invariant 
of the second degree f2l| , p2| . The two knot invariants are practically useful for computer 
simulation of random polygons with a large number of polygonal nodes. In fact, it has been 
demonstrated in Ref. |H that the Vassiliev invariant V2{K) can be calculated not only in 



polynomial time but also without using large memory area. 
B. Characteristic length of random knotting N c (r) 

For a given knot K, we consider the probability Pk{N,t) that the topology of an N- 
noded self-avoiding polygon with cylinder radius r is given by the knot type K. We call 
it the knotting probability of the knot K. Let us assume that we have Mk self-avoiding 
polygons with a given knot type K among M samples of cylindrical SAPs with radius r. 
Then, we evaluate the knotting probability Px(N,r) by Px(N,r) = Mk/M. 

For the trivial knot, the knotting probability P tr i V (N,r) for the cylindrical SAPs is given 

by 

P triv (N, r) = C triv exp(-N/N c (r)) . (1) 



Here the estimate of the constant Ctriv is close to 1.0 ||13|| . We call N c (r) the characteristic 
length of random knotting. It is also shown in Ref. that N c (r) can be approximated by 
an exponential function of r: 

iV c = iV c (0)exp(7r). (2) 
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The best estimates of the two parameters N c (0) and 7 are given by N c (0) = 292 ± 5 and 



7 = 43.5 ±0.6 13 



For several knots, it is shown [|14| that the knotting probability Pk(N, r) of a knot K is 
given by 

\ m ( K ) 



P K (N,r) = C K 



N K (r) 



exp(-N/N K (r)) 



(3) 



It is numerically suggested in Ref. JL4J] that Nx(r) should be independent of K: Nk(t) 
N c (r), and also that the constant Ck should be independent of the cylinder radius r. 



C. Mean-squared gyration radius with a topological constraint 

The mean square radius of gyration R 2 of a self-avoiding polygon is defined by 

1 - 

R2 = ^-2 E <0Rn-i? m f> • (4) 
n,m=l 

Here R n is the position vector of the nth segment (or the nth node) and < • > denotes the 
ensemble average, which is taken over all possible configurations of the self-avoiding polygon. 

Suppose that we have M self-avoiding polygons. Then, we evaluate the mean square 
radius of gyration R 2 by the sum: R 2 = Y^iLi R-i/M, where R 2 denotes the gyration radius 
of the ith SAPs in the given M SAPs. 

Let us define the mean square radius of gyration R 2 K for such self-avoiding polygons that 
have a given knot type K: 

1 M K 

n 8=1 

where R 2 Ki denotes the gyration radius of the zth self-avoiding polygon that has the knot 
type K. In terms of R 2 K , R 2 is given by R 2 = Y,k RrM k /M. 

In Fig. 1, the estimates of the mean square radius of gyration R 2 are plotted against the 
number A^ of nodes in a double-logarithmic scale, for the cylindrical SAPs with r = 0.003 
and r = 0.03. We may confirm the standard asymptotic behaviors of the mean-squared 
gyration radius R 2 in Fig. 1. Here we remark on an effective exponent zv eff , which is defined 
through the power- law approximation: R ~ N UcB . It is shown in [^] that the estimate of 
the effective exponent i^af for the cylindrical SAPs with radius r is consistent with that of 
the cylindrical SAWs with radius r. 
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III. FINITE-SIZE BEHAVIORS OF R 2 FOR SOME KNOTS 



Let us discuss simulation results on the mean square radius of gyration R\ for the 
cylindrical SAPs with a knot K and of radius r. For two prime knots (the trivial and 
trefoil knots) and a composite knot (the double-trefoil knot, 3ijj3i), we have investigated 
the mean-squared gyration radius R 2 K under the topological constraint in the range of the 
number N satisfying 21 < N < 1001, and for 14 different values of cylinder radius r. 

The gyration radius R 2 K can approximately given by some power of N . In Fig. 2, double- 
logarithmic plots of R 2 K versus N are given for the trivial and trefoil knots, with two values 
of cylinder radius: r = 0.003 and 0.03. We see that all the double-logarithmic plots of 
Fig. 2 fit to some straight lines. We note that for other values of cylinder radius r, several 
double- logarithmic plots of R 2 K versus N are explicitly shown in Ref. p3| . 

With the number N fixed, R\ should increase with respect to the radius r for any knot. 
In Fig. 2, closed squares for r = 0.03 are located higher in the vertical direction than closed 
circles for r = 0.003, through the whole range of N both for the trivial and trefoil knots. 

A. Ratio R 2 K /R 2 and the effective expansion under the topological constraint 

Let us now consider the ratio of R\ to R 2 for a given knot K. If the ratio is larger (smaller) 
than 1.0, then the average size of SAPs with the knot K is relatively larger (smaller) than 
that of no topological constraint. We say that the SAPs with the knot K is effectively 
more (less) expanded. In Fig. 3, the ratio R 2 K /R 2 versus the number N is plotted in a 
double-logarithmic scale for the trivial and trefoil knots. Here, we have depicted only the 
case of r = 0.003 among many sets of the cylindrical SAPs with the 14 different values of 
cylindrical radius. 

For the trivial knot, we see in Fig. 3 that the ratio R 2 riv /R 2 is greater than 1.0 when 
N > 50. Thus, the average size of the ring polymers with the trivial knot enlarges under 
the topological constraint. It gives a typical example of effective expansion. 

In Fig. 3, the graph of the trivial knot is convex downwards: the ratio R 2 riv /R 2 is almost 
constant with respect to N for small N such as N < 100; for N > 300 the ratio R 2 riv /R 2 
increases with respect to N with a larger gradient, and the graph can be approximated by 
a power law such as R\ riv jR 2 oc N v <& ". Here the symbol i/^ v denotes the effective exponent 
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for the trivial knot. We note that the characteristic length N c (r) is approximately given by 
300 for r = 0.003. Thus, we may say that the power law behavior is valid for N > N c (r). 

For the trefoil knot, the graph can be approximated by a power of N such as R 2 re oc N v >& 
through the range of 100 < iV < 1001. Here the symbol u^f denotes the effective exponent 
of the trefoil knot. In Fig. 3, we find that when N < 100 the ratio R 2 re /R 2 is smaller 
than 1.0, while it is larger than 1.0 when N > 300. Thus, when N is small, the topological 
constraint of the trefoil knot gives effective shrinking to ring polymers, while it does not 
when N is large. For a nontrivial knot K, we expect that the ratio R 2 K /R 2 is less than 1.0 
when N is small, while it can be larger than 1.0 when N is large. 

The properties of the ratio R 2 K /R 2 discussed in the last three paragraphs are consistent 
with the simulation results of Gaussian random polygons |24| . We have found for the random 
polygons that the double- logarithmic graph of R\jR 2 versus iV is given by a downward 
convex curve for the trivial knot, while it is given by a straight line for the trefoil knot 
and also for other several nontrivial knots; for the nontrivial knots investigated, the ratio 
R\ I R 2 is given by some power of iV such as N v <& . Thus, there are indeed many important 
properties valid both for the simulation of the Gaussian random polygons and that of the 
cylindrical SAPs with a small radius such as r = 0.003. 

The observations derived from Fig. 3 should be valid particularly for finite-size systems. 
Admitting that N is finite, we can only understand that the Gaussian random polygons 
and the cylindrical SAPs have the similar topological properties in common. If we discuss 
asymptotic behaviors, SAPs and random polygons should be quite different. However, if we 
consider such properties that are valid for finite N, then they can hold both for SAPs with 
small excluded volume and random polygons that have no excluded volume. 

Let us discuss again the convexity of the graph of the trivial knot, which has been observed 
in Fig. 3. We consider how the convexity depends on the radius r. In Fig. 4, the graphs of 
the ratio R 2 riv /R 2 versus N are given in a double- logarithmic scale for four different values 
of cylinder radius r. Then, we see that the graph with r = 0.05 is less convex than that of 
r = 0.003. Thus, the convexity in the graphs of the effective expansion for the trivial knot 
should be valid only when cylinder radius r is small. 

Let us assume that the convexity of the graphs of R 2 riv /R 2 for the small r case should 
correspond to a crossover behavior of R\ riv jR 2 with respect to N. Then, the crossover 
behavior could be related to that of Gaussian random knots, which is recently discussed by 
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Grosberg |18| for Gaussian random polygons. We can discuss the convexity of the double- 
logarithmic graph of R 2 riv /R 2 versus N, taking an analogy with the crossover of the Gaussian 
random knots. Thus, we call the convexity of the trivial knot in Fig. 3 the crossover, 
hereafter in the paper. 

For the non-trivial knots investigated, we do not see any crossover in the graph of R 2 K / R 2 
versus N. For instance, for the 4i and 3ij}3i knots, the slope of the graph near iV ~ N c (r) 
is straight in the double-logarithmic scale. The crossover at iV ~ N c (r) should be valid only 
for the trivial knot. 



B. The plateau in the graph of R 2 K /R 2 versus N for large N 

We discuss how the ratio R 2 K / R 2 depends on the number N, considering both the excluded 
volume effect and the effective expansion due to the topological constraint. In Fig. 5, the 
graphs of the ratio R\jR 2 versus N for different values of cylinder radius r are shown in 
linear scales: (a) for the trivial knot; (b) for the trefoil knot. 

Let us first consider the large N behaviors of the graphs shown in Fig. 5 for the trivial and 
trefoil knots. The graphs of R 2 K /R 2 versus iV have a common tendency that they become 
constant with respect to N when iV is very large. It is particularly the case for the larger 
values of cylinder radius r such as r = 0.03 and 0.05. They approach horizontal lines at 
some large values of N. When r is small such as r = 0.003, the graph becomes flat only for 
large N, as shown in Fig. 5. 

From the flatness of the graphs of R 2 K /R 2 for large N, it follows that the power law 
behavior: R\ r jR 2 oc iV* does not hold when N is very large. In Fig. 3, we have discussed 
that the ratio R 2 re /R 2 versus the number iV can be approximated by the power law for 
r = 0.003 through the range of 100 < N < 1001. However, the power-law approximation 
should be valid only within some finite range of N. 

Let us discuss other finite- N behaviors of the ratio R 2 K jR 2 . For the trefoil knot, the ratio 
R 2 re /R 2 is less than 1.0 when N is small; it approaches or becomes larger than 1.0 when N 
is large enough. When cylinder radius r is small such as r = 0.003 and r = 0.01, the ratio 
R 2 re /R 2 is clearly greater than 1.0 when N is large enough. When r is small, there should 
be a critical value N critica i such that R 2 re /R 2 < 1.0 for N < N criticah and R 2 re /R 2 > 1.0 
for iV > N critical- Furthermore, we have a conjecture that the critical value N critica i should 



S 



be roughly equal to the characteristic length N c (r) of random knotting. It seems that the 
conjecture is consistent with the graphs of Fig. 5 (b). 

Let us discuss the conjecture on N critica i, explicitly. In Fig. 5 (b), we see that for 
r = 0.003, the ratio R 2 tr jR 2 becomes 1.0 roughly at N = 300, and also that for r = 0.01, 
the ratio R 2 r jR 2 is close to 1.0 roughly at N = 400. The observations are consistent 
with the estimates of N c (r) in Ref. @: N c (r) = (2.72 ± 0.06) x 10 2 for r = 0.0 and 
N c {r) = (4.72 ± 0.14) x 10 2 for r = 0.01. Thus, the consistency supports the conjecture on 

^critical • 



C. Decrease of the topological effect under the increase of the excluded volume 

The effect of a topological constraint on the gyration radius decreases when the excluded 
volume increases. There are two examples: the decrease of ratio R 2 K /R 2 with respect to 
cylinder radius r while N being fixed, and the disappearance of the crossover for the trivial 
knot shown in Figs. 3 and 4. 

Let us first discuss how the excluded-volume can modify the effective expansion due to 
the topological constraint. As we clearly see in Fig. 5, the ratio R 2 K jR 2 decreases as cylinder 
radius r increases with iV fixed, both for the trivial and trefoil knots. Thus, the effective 
expansion of SAPs under the topological constraint becomes smaller when the excluded 
volume becomes larger. 

It is quite nontrivial that the effective expansion given by the ratio R 2 K /R 2 decreases as 
cylinder radius r increases. In fact, the value of R 2 K itself increases with respect to r, as 
we have observed in Fig. 2. Furthermore, one might expect that the effective expansion 
due to a topological constraint should also increase with respect to cylinder radius r, simply 
because the average size of ring polymers with larger excluded volume becomes larger, as 
observed in Fig. 1. However, it is not the case for the ratio R 2 K /R 2 . 

Let us now discuss the crossover behavior of the trivial knot again, from the viewpoint 
of the competition between the topological effect and the excluded volume effect. Here we 
recall that the crossover has been discussed in §3. A with Figs. 3 and 4. Here we regard the 
crossover as a characteristic behavior derived from the topological constraint of being the 
trivial knot. 

As a working hypothesis, let us assume that the crossover should occur at around the 
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characteristic length N c (r). Recall that N c (r) is larger than 1000 for r = 0.03 and 0.05, as 
we have estimated: N c (r) « 1200 for r = 0.03, and N c (r) « 2600 for r = 0.05. If the above 
hypothesis would be valid, then the graphs for r = 0.03 and 0.05 should also be convex. In 
Fig. 4, however, we see no change in the gradient of the graph of R 2 riv /R 2 versus iV for 
r = 0.03 or 0.05. The assumed crossover of the trivial knot does not appear for r = 0.03 or 
0.05. We may thus consider that the crossover as a topological effect is diminished by the 
excluded volume effect when r > 0.03. 

D. Characteristic length of random knotting N c (r) and the effective expansion 

In terms of the characteristic length N c (r), we can explain some properties of the effective 
expansion of cylindrical SAPs under a topological constraint. Here we recall that the ratio 
R\j R 2 describes the degree of the effective expansion under the topological constraint of a 
knot K. 

We first consider the case when the characteristic length N c (r) is very large. Let us show 
that the ratio R 2 riv /R 2 should be close to 1.0 for JV •< N c (r). First, we recall that the 
probability P triv (N) of the trivial knot decays exponentially with respect to the number N 
of polygonal nodes: P triv (N) = exp(— N / N c (r)) . If N/N c (r) is very small, the probability 
P triv (N) is close to 1.0, i.e., almost all SAPs have the trivial knot. Then, the mean-squared 
gyration radius with no topological constraint R 2 should be almost equal to that of the 
trivial knot R 2 riv . Consequently, the ratio R 2 riv /R 2 should be close to 1.0. 

When r > 0.05, the characteristic length N c (r) is larger than 2600. Then, the trivial 
knot is dominant among the possible knots generated in SAPs with N < 1000. Thus, R 2 riv 
should almost agree with R 2 , which is the mean-squared gyration radius of SAPs under no 
topological constraint. There is no effective expansion under the topological constraint: the 
R 2 riv /R 2 is close to 1.0. 

Let us next consider the case when the characteristic length N c (r) is small or not large. 
Then we show that the mean square radius of gyration of SAPs with the trivial knot R 2 triv 
should be larger than that of no topological constraint R 2 for N > N c (r). In fact, various 
types of knots can appear in a given set of randomly generated SAPs of the cylinder radius 
r, since the probability of the trivial knot P triv (N) is exponentially small for iV > N c (r) . 
We note that the fraction of nontrivial knots is given by 1 — exp(—N/N c (r)). Thus, it is not 
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certain whether the ratio R 2 riv /R 2 is close to the value 1.0 or not. However, we may expect 
that the ratio R 2 riv /R 2 should be indeed larger than 1.0. Here we consider the following 
points: when N > N c (r), the majority of SAPs generated randomly should have much more 
complex knots than the trivial knot; the mean square radius of gyration of iV-noded SAPs 
with a very complex knot should be much smaller than that of the trivial knot. 

The explanation on the effective expansion discussed in the above is completely consistent 
with the simulation results, as having been discussed in §3, in particular, through Figs. 3, 
4 and 5. 

IV. ASYMPTOTIC BEHAVIORS OF R 2 K 
A. The exponent of R 2 K 

Let us discuss an asymptotic expansion for the mean square radius of gyration of cylin- 
drical SAPs with a given knot K. Here we assume that R 2 K can be expanded in terms of 
1/N consistently with renormalization group arguments. Then, the large N dependence of 
R 2 K is given by 



Here, the exponent vk should be given by that of self-avoiding walks: vk = ^saw- in order 
to analyze the numerical data systematically, however, we have introduced Vk as a fitting 
parameter. Thus, for the ratio R 2 K /R 2 , we have the following expansion: 



Here we have put Avr fitting parameter. 

We have analyzed the data for the three different knots: the trivial, trefoil and 3i(j3i 
knots, applying the expansion (0) to the numerical data of R 2 K /R 2 for N > 300 . The best 
estimates of the three parameters are given in Tables 1, 2 and 3 for the trivial, trefoil and 
3i(j3i knots, respectively. 

Let us discuss the best estimates of the difference of the exponents: Avk- We see in 
Tables 1 , 2 and 3 that all the results of Av^ suggest that they should be given by 0.0, with 
respect to the confidence interval. Let us examine the best estimates more precisely. It is 



R\ = A K N 2VK 1 + B K N- A + 0(1/N) . 



(6) 



R\jR 2 = (A k /A)N 2Auk 1 + (B K - B)N~ A + OQ./N) 



(7) 
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rather clear from Tables 1, 2 and 3 that for a given cylinder radius r, the best estimates of 
Auk are independent of the knot type. 

There is another evidence supporting that Ai>k = 0.0 for the trivial and trefoil knots. 
Let us consider the plots of the ratio R 2 K /R 2 versus N in Fig. 5 for the trivial and trefoil 
knots. We recall that the graphs are likely to approach some horizontal lines at some large 
N. The tendency of the graphs becoming flat for large N suggests that R K and R 2 should 
have the same exponent, i.e., vsaw- 

From the two observations, we conclude that the difference of the exponents is given by 
0.0: Avk = 0.0 for any value of r. There is thus no topological effect on the scaling exponent 
defined in the asymptotic expansion of R 2 K . 

B. Amplitude ratio Ak/A 

Let us now consider the amplitude Ak of the asymptotic expansion (0). In Tables 1, 2 
and 3, the best estimates of the ratio Ak/A are larger than 1.0 for the three knots, when r 
is small. The observation must be important. In fact, if the amplitude ratio Ak/A is larger 
than 1.0 in the asymptotic expansion (|7|), then R\ is larger than R 2 for any large value 
of N. It might seem that the consequence is against the standard thermodynamic limit of 
statistical mechanics. 

However, there is a clear evidence for the observation that Ak/A > 1.0 for some small 
values of cylinder radius r. In fact, the graphs of the ratio R 2 K /R 2 versus N are monotonically 
increasing with respect to N, as we see in Figs. 3, 4 and 5. It is clear that the graphs with 
the smaller values of cylinder radius r are larger than 1.0 when N is large. This observations 
of Figs. 3, 4 and 5 confirm that Ak/A > 1.0 when cylinder radius r is small. Thus, we may 
conclude that the topological constraint gives an effective expansion also to asymptotically 
large cylindrical SAPs when the radius r is small. 

The value of Ak/A decreases with respect to the radius r for the three knots. We see 
it in Tables 1 to 3, where the best estimates of Ak/A are listed. It is also consistent with 
the fact that the ratio R 2 K /R 2 decreases with respect to r, which we have discussed in §3.C. 
However, the decrease of Ak/A is quite nontrivial, since the mean-squared gyration radius 
R\ itself increases with respect to r, for the trivial and trefoil knots, as shown in Fig. 2. 
Here we recall in Fig. 1 that the gyration radius under no topological constraint R 2 increases 
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with respect to r. 

From the viewpoint of asymptotic behaviors, we have shown that the effective expansion 
derived from the topological repulsion decreases with respect to cylinder radius r. We have 
also discussed that R\ is larger than R 2 for any large value of N, when cylinder radius r is 
small. 

C. The r-dependence of the amplitude ratio 

Let us discuss the r dependence of the amplitude ratio Ak/A, more quantitatively. For 
this purpose, we analyze the data of R\/R 2 versus N again, assuming vk = v in eq. (0). 
We evaluate the amplitude ratio A K /A by the following formula: 

R 2 K /R 2 = a K (l + (3 K N- A + 0{l/N)) . (8) 

Here we have replaced with ax and /3k, Ak/A and Bk — B in (0), respectively. Here we 
have also introduced a technical assumption: A = A^=0.5 in (|7|). 

We have obtained the numerical estimates of ax, applying the fitting formula (§) to the 
data of R 2 K /R 2 with N > 300 . The estimates of a.K versus r are shown in Fig. 6 in the 
double-logarithmic scale for the trivial, trefoil, and 3ijj3i knots. To be precise, the values of 
ax are a little larger than those of Ak/A given in Tables 1 , 2 and 3. 

The estimate of the parameter becomes close to the value 1.0 when cylinder radius r 
is large enough. Furthermore, it is suggested from Fig. 6 that should be independent of 
the knot type. In fact, the data points for the trivial, trefoil and the double-trefoil (3i#3i) 
knots overlap each other. These two observations are consistent with the simulation result 
of the self-avoiding polygons on the lattice JF|, |TTJ . 

Interestingly, we see in Fig. 6 that the ratio o:k decreases monotonically with respect to 
the cylinder radius r. For the data with 0.001 < r < 0.01, we find that ctK is roughly ap- 
proximated by a decreasing function of r such as ok = aor^ exp(— tpr), with cto=1.00±0.12, 
0=-O.O5±O.O2 and ^=5.78±4.79. The x 2 value is given by 1. 

V. DISCUSSION 

With some graphs in the N — r plane, we can illustrate the finite-size behaviors of the 
ratio R 2 K /R 2 discussed in §3. We recall that the topological effect has played a central role 
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as well as the excluded- volume effect. Thus, we consider two lengths with respect to the 
number iV of polygonal nodes: the characteristic length of random knotting N c (r) and the 
"excluded-volume length" N ex (r). When iV > N ex (r), the excluded-volume effect should be 
important to any iV-noded SAP with radius r . 

We define N ex (r) by N ex (r) = 1/r 2 . The derivation is given in the following. We first 
note that the parameter z of the excluded-volume is given by z = Const. \fNB /£ 3 oc N 1 / 2 r, 
where the cylindrical segments have the diameter d and the length £, and the second virial 
coefficient B of a polymer chain is given by £ 2 d [^5[. Here we also note that the ratio d/£ 
corresponds to the radius r of the cylindrical SAPs. We may consider that when z « 1, the 
excluded volume can not be neglected. Thus we have the number N ex (r) from the condition: 
^N ex (r)r = 1. 

We consider two graphical lines in the N — r plane: A" = N ex (r) and A" = N c (r). In Fig. 
7, the vertical line of the diagram expresses the r-axis and the horizontal one the A^-axis. 
The graph N c (r) = N reaches the A" axis at A^ = A^ c (0) ~ 300. Here we recall that the 
function N c (r) is given by eq. (fj): N c (r) = N c (0) exp (77*). There is a crossing point for the 
two curved lines. The coordinates of the crossing point is approximately given by N* = 1300 
and r* = 0.03. For a given simulation of the ratio R 2 K /R 2 with a fixed radius r, we have a 
series of data points located on a straight line parallel to the A" axis. 

Let us first consider the case of small values of r such as r = 0.003 and r = 0.01. From the 
simulation of §3, it is shown that the effective expansion due to the topological constraint 
is large. This is consistent with the following interpretation of the N — r diagram: if we 
start from the region near the r axis and move in the direction of the A" axis, then we 
cross the line A" = A r c (r) before reaching another one viV r = 1; thus, we expect that the 
excluded-volume remains small when the topological effect becomes significant. 

The above explanation should be consistent with the observation that the crossover of 
the trivial knot occurs near A" = N c (r) for small values of r. Here we recall Figs. 3 and 4. 
When r is very small, then we cross the line of N c (r) = N almost at N c (0) ~ 300. 

When radius r is large such as r = 0.03 and 0.05, it is shown in §3 through simulation 
that the effective expansion is small: the ratio R 2 K / R 2 is close to 1.0. In the N — r diagram, 
when we move rightwards from the region near the r axis with r fixed, we cross the line 
\/N r = 1 before reaching another line N c (r) = N. Thus, the effective expansion as the 
topological effect should be small. 
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Finally, we should remark that some important properties of R 2 K of cylindrical SAPs with 
radius r have been discussed systematically through scaling arguments with the blob picture 
by Grosberg |[26|| . In the note p6 |, the characteristic length N c (r) and the excluded- volume 
parameter z are explicitly discussed in the N — r diagrams. It would thus be an interesting 
future problem to investigate how far the predicted properties of R 2 K are consistent with 
simulation results. 
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FIG. 1: Double-logarithmic plots of the mean square radius of gyration under no topological 
constraint R 2 for cylindrical SAPs versus the number N of polygonal nodes for radius r = 0.003 
and 0.03 depicted by closed circles and squares, respectively. 

FIG. 2: Double-logarithmic plot of R 2 K versus N with r = 0.003 and 0.03 shown by closed circles 
and squares, respectively: (a) for the trivial knot; (b) for the trefoil knot. 

FIG. 4: Double-logarithmic plots of the ratio R 2 riv /R 2 versus N for r =0.003, 0.01, 0.03 and 0.05 
shown by closed circles, squares, diamonds and triangles, respectively. 



FIG. 5: Graphs of the ratio R 2 K /R 2 versus the number iV in linear scales for r = 0.003, 0.01, 0.03 
and 0.05 shown by closed circles, squares, diamonds and triangles, respectively: (a) for the trivial 
knot and (b) for the trefoil knot. The same data points are shown in both Fig. 4 and Fig. 5 (a). 



FIG. 6: Double-logarithmic plots of the amplitude ratio ax versus cylinder radius r for the trivial, 
trefoil and double-trefoil (3i#3i) knots shown by closed circles, squares and triangles, respectively. 
For the double-trefoil knot, the data points for 0.001 < r < 0.01 are shown. 



FIG. 7: N — r diagram. Graphs of N = N c (r) and N = N ex {r) are shown by two curved lines. 
The arrows (a) and (b) suggest the series of the data points of Fig. 5 for r = 0.01 and r = 0.005, 
respectively. All the data points in the paper are located in the shaded area. 

FIG. 3: Double-logarithmic plots of the ratio R 2 K /R 2 versus N for cylindrical SAPs with r = 0.003. 
R 2 riv /R 2 and R 2 re /R 2 are shown by closed circles squares, respectively. 
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TABLE I: Fitting 


; parameters Ak/A, Bk 


— B and Auk 


versus cylinder radius r: 


for the trivial 


knot. 










r 


A k /A 


Bk-B 


2Av K 


x 2 


0.001 


1.313±1.285 


-2.587±4.211 


0.003±0.123 


18 


0.002 


1.235±1.152 


-2.389±4.117 


0.009±0.116 


12 


0.003 


1.213±1.065 


-2.317±3.912 


0.009±0.109 


3 


0.004 


1.228±1.062 


-1.982±3.973 


0.003±0.107 


3 


0.005 


1.170±0.983 


-1.684±3.985 


0.007±0.104 


4 


0.006 


1.207±0.921 


-2.204±3.464 


0.005±0.095 


4 


0.007 


1.159±0.891 


-1.633±3.684 


0.005±0.095 


3 


0.01 


1.106±0.836 


-1.090±3.809 


0.005±0.092 


3 


0.02 


1.065±0.686 


-0.699±3.384 


0.003±0.078 


1 


0.03 


1.063±0.628 


-0.518±3.166 


-0.001±0.071 


2 


0.04 


1.043±0.590 


-0.353±3.076 


-0.001±0.068 


1 


0.05 


1.010±0.554 


-0.143±3.039 


0.002±0.066 


1 


0.06 


1.020±0.551 


-0.103±2.997 


-0.001±0.065 


1 


0.07 


1.013±0.531 


-0.187±2.898 


-0.001±0.060 


1 
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TABLE II: Fittinj 


I parameters Ak/A, Bk 


— B and Auk 


versus cylinder radius r: 


for the trefoil 


knot 










r 


A k /A 


Bk-B 


2Avk 


x 2 


0.001 


1.286±0.970 


-4.440±2.784 


0.014±0.096 


3 


0.002 


1.215±0.918 


-4.093±2.906 


0.015±0.095 


10 


0.003 


1.202±0.905 


-3.562±3.054 


0.011±0.094 


11 


0.004 


1.176±0.872 


-3.423±3.069 


0.012±0.093 


16 


0.005 


1.176±0.845 


-3.461±2.964 


0.010±0.090 


6 


0.006 


1.113±0.807 


-3.174±3.091 


0.015±0.090 


7 


0.007 


1.084±0.765 


-3.219±2.996 


0.019±0.088 


3 


0.01 


1.103±0.743 


-3.220±2.870 


0.013±0.084 


1 


0.02 


1.068±0.765 


-2.326±3.353 


0.005±0.088 


2 


0.03 


1.058±0.790 


-2.262±3.531 


0.003±0.091 


4 


0.04 


1.003±0.835 


-2.043±4.034 


0.007±0.101 


3 


0.05 


1.007±0.883 


-2.422±4.119 


0.007±0.107 


4 


0.06 


1.029±0.975 


-2.900±4.274 


0.005±0.116 


3 


0.07 


0.998±1.197 


-1.923±5.915 


0.002±0.146 


2 



TABLE III: Fitting parameters Ak/A, Bk — B and Avk versus cylinder radius r: for the double- 
trefoil knot (3ijJ3i). 



r 


Ak/A 


Bk - B 


2Av K 


x 2 


0.001 


1.269±1.158 


-5.203±3.255 


0.012±0.116 


7 


0.002 


1.224±1.077 


-5.203±3.113 


0.016±0.112 


7 


0.003 


1.158±1.090 


-4.371±3.662 


0.016±0.118 


1 


0.004 


1.149±1.054 


-4.866±3.401 


0.018±0.116 


9 


0.005 


1.137±1.008 


-4.851±3.297 


0.016±0.112 


4 


0.006 


1.096±1.002 


-4.745±3.476 


0.021±0.115 


1 


0.007 


1.061±1.043 


-3.819±4.091 


0.020±0.122 


3 


0.01 


1.076±1.369 


-3.563±5.287 


0.012±0.159 


4 
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